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1 Introduction 

One of the major purposes of control theory is the analysis and the design of 
feedback control systems. A great number of both analytical and computational 
methods have also been adopted for continuous and discrete systems, [5j,|14j. 
In this paper we examine the feedback stabilization problem for a large variety 
of continuous nonlinear systems. Our methodology resembles that of the applied 
computational tools used to solve control problems jTHj. Specifically, we consider 
a system of the form x = $(x, u), (E), where the components of $i, are mul- 
tivariable polynomials. According to the theory, if the linearization of (S) at the 
origin is asymptotically controllable (i.e., all the uncontrollable eigenvalues have 
negative real parts), then the nonlinear system is locally asymptotically stabiliz- 
able This means that we can find a feedback-law that makes the closed-loop 
system asymptotically stable to the origin. The aim of this paper is to calculate 
those feedback-laws computationally by means of certain symbolic algorithms. 
Concretely, we seek for nonlinear state-feedbacks of the form u = a(x), where 
a(x) consists of multivariable polynomials, so that a proper Lyapunov function 
becomes negative. To solve the problem we develop the following algorithms: 
The Formal Algorithm. This algorithm allows us to write a polynomial p as 
follows: 

P = ci(H^.,<.,^)xf ^ ■ [1^2,1,13:1 + x^]^^'^ ■ [M/3,1,1^1 + M^3,2,ia;2 + x^]^''' ■ ■ ■ 

■ ■ ■ [Wn,l,lXi + WnXlX2 + " " " + + 
+C2{W,,.,^)x{''' ■ [W2,l,2Xl + X^Y^'^ ■ [1^3,1,2X1 + ^3,2,20:2 + X^^^^^ ■ ■ ■ 
■ ■ ■ [Wn,l,2Xl + I^n,2,2a;2 + " " " + X^Y^'^ + ■■■ 
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■ ■ ■ [WnXkXl + Wn,2,kX2 + " " " + + i?w(Xi) (1) 

where the exponents ia,h are specific positive whole numbers, the quantities Wi^a,<p 
are undetermined parameters that can take certain values, Cj(Wi,o-,(/.) the coeffi- 
cients depending on the parameters Wi^a,<f> and i?w(xi) a polynomial of the single 
variable xi, called the remainder. Equation (0) is called the Formal-Linear-Like- 
Factorization of p and appeared firstly in [15j. The essential tool of this methodol- 
ogy is a continuous reduction of p, by means of an Euclidean division. This makes 
our method similar to others [HI, [H]. Yet in our methodology we only deal with a 
specific polynomial, in addition to the Grobner basis that works with a polynomial 
ideal. 

The POS- Algorithm. This algorithm checks the positivity of a multivariable 
polynomial p. This task can be achieved by giving to Wi^a,4> values such that all the 
non-square terms to be eliminated and the square terms have positive coefficients. 
A lot of work has been done in this direction and some of the results can be 
found in the following PP,[Sl,J2l! to mention but a few. In recent years there 
has been a strong algorithmic approach to the whole issue, as can be seen, for 
instance, in jl],|n],jH]- A number of those algorithms dealt with the problem of the 
existence of a sum of squares decomposition |2] , [Z] , IIHI • That is when we are able 
to write a polynomial p as p = fii^i, X2, ■ ■ ■ , a;„). In our sum of squares, p = 
X]j fii^ii X2, ■ ■ . , Xn), the quantities fi{xi,X2, ■ ■ ■ , Xn) are " linear " polynomials of 
the variables Xi,X2, ■ ■ ■ , Xn, while in other approaches [TU], fi{xi, X2, ■ ■ ■ , Xn) can 
be arbitrary. Moreover, our methodology provides not a specific sum of squares 
but a whole class, depending on the various proper values of the parameters Wi^a,-^- 
The Feedback-GAS- Algorithm. This algorithm accepts as input the poly- 
nomials <I>j and a Lyapunov Function L, and calculates the feedback connection 
u = a(x). To achieve this we assume that a(x) consist of a multivariable poly- 
nomial with parametric coefficients (we denote them by A^^^^^ . J. Our aim is 

to determine those values of the parameters Wi^„^^ and A^-'^\^ that make the 
Lyapunov function negative along the trajectories of the closed-loop system. This 
guarantees the stability of the origin. The main merits of our method are the 
following: 

1) It constitutes a pragmatic computational method. Indeed, the feedback-laws 
are derived from symbolic computational algorithms. Appropriate software has 
been created for this purpose, and all the examples presented in the current paper 
have been studied with the aid of this software. 

2) It provides us with a whole class of admissible controllers, as opposed to only a 
single one. 

3) It works even when the linearization of the nonlinear system has an uncontrol- 
lable eigenvalue with zero real part, |14j . 
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4) Since the Lyapunov function, upon construction, is negative for every point, 
the stabihty is global. 

The paper is divided into two parts. The first part is devoted to the presentation 
of the algorithms used throughout the paper. The second part is devoted to the 
specific control application. R and Ti^ will denote the sets of real numbers and 
positive integers correspondingly. 

2 The Algebraic Background 

Let R be the set of real numbers and Xi, X2, . . . , x„ n- variables. An expression 
of the form p = X]a=i c\x°^'^x^2 '^ ' ' "^^n"'^, where ca G R and some of the expo- 
nents e Z+ are not equal to zero, is called a polynomial in xi, X2, . . . , Xn with 
real coefficients or, for short, a real polynomial. The set of all real polynomi- 
als in xi, X2, ... , Xn is denoted by R[xi, X2, . . . , xj. An element x^^'^Xa^'^ ■ ■ ■ Xn"'^ 
is called a monomial and an element called a term. Let 

0n,A = x"^'^X2^'^ ■ ■ ■ x^n'^ and = x"^'''x2^''' ■ ■ ■ Xm"'** be two mouomials. We de- 
fined the lexicographical order among monomials jD], as follows: we say that 0„,a 
is ordered less than (j)m,^l, denoted by 0n,A -< 4>m,n, if either n < m or n = m and 
o-n,x < cim,^i- In other words, the monomials are ordered as follows: Xi ^ ■ ■ ■ -< xj" -< 
• • ■ ^ X1X2 -< • • ■ ^ X1X2 -<■■■-< X1X2X3 -< ■ ■ ■ Let p be a given polynomial. Or- 
dered lexicographically, the term that corresponds to the maximum monomial is 
called the maximum term denoted by maxterm{p) . Throughout the paper, a vari- 
able taking values in R is called an undetermined parameter. The set of 
undetermined parameters is denoted by W. 

Let p G R[xi,X2, . . . ,Xn] be a polynomial with n-variables and W = 
a set of undetermined parameters, taking values in R. A Formal-Linear-Like- 
Factorization of p is an expression of the form: 

k 

P=Y^ C,^{Wi,^M''" ■ [W2,l,t.Xi + X2y''^ ■ [W3,l,^,Xl + W3,2,^X2 + Xg]^'^'^ • • • 
■ ■ ■ [Wn,l,^,Xl + Wn,2,^X2 + " " " + Xn]^"'" + Rw{xi) 

where the coefficients Cf^iWi^a,-^) are polynomial functions of the parameters Wi^a^^ 
and the remainder Ry^{xi) is a polynomial only of the single variable Xi, with 
coefficients depending on the parameters Wi^^^^, too. Some of the exponents 
ji,M)i2,/i) • • • )in,/i G Z"*", may be equal to zero. The Formal- Linear- Like- Factorization 
of p is denoted by FormalLF[p]. In a number of instances, a Formal-Linear-Like- 
Factorization of p, is written for short as FormalLF[p] = X]^=i '^f^^i^,'/^ ' ^'2,1^ ' 

^3,'/^ ■ ■ ■ LnJl^ + -R, where L„^^ = W^^i^^xi + W„^24iX2 H h x^, a = 1, . . . , n and 

c^,i? are the abbreviations of c^{Wi^a,^p) and i?>v(xi). 
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Example 2.1 We have the polynomial p — 5xi — 7xiX2 + llxix^. The Formal- 
Linear-Like-Factorization of p is, p = llxi(H^3^i^iXi+ 1,23^2 + 3:3)+ (—7 — 

11^^3,1,2) Xi{W2,l,2Xl + X2)+ 5X1 + (7H/2,1,2- 11^^3,1,1+ llVI/2,l,2W/'3,2,l)a;? • 

The following theorem deals with the question of uniqueness of the Formal- Linear- 
Like- Factorization. 

Theorem 2.1 For a given polynomial p e R[xi,X2, . . . ,Xn]; the FormalLF[p\ is 
unique, under the assumption that the parameters W — {Wi j,p} are considered as 
constants. 

Proof: If p has only xi-terms, the proof of the theorem is trivial, with = 
and R — p. Let us suppose that p has at least one term other than the Xi-terms. 
Let us further suppose that Ax'J^'''x2^''' ■ ■ ■Xn" '' is its maximum term. This term 
appears also in the product: ChLl^'/^ L'l^'j^ ■ ■ ■ L'^^h ■ equating their coefficients, 
we calculate the quantity c/j, uniquely, actually c/j = A. Repeating the same 
procedure for the term with the next higher order, we find an expression for the 
"next" coefficient c/j_i. Since this expression is a function of Ch and some of the 
parameters Wi^cr,ip are considered as constants, we conclude that Ch-i is also defined 
uniquely. By induction, we finally get coefficients c^, ^ = 1, . . . , A;, all of which are 
uniquely determined. The polynomial R consists only of Xi-terms. These terms 
arise either from the polynomial p or from the products c^L^i^'l^ ■ L^^'l^ ■ Lf '^^ ■ ■ ■ Ln,'^ . 
The unique determination of the coefficients entails the uniqueness of i?, and 
the theorem has been proved. □ 

What is of interest is the issue of constructing the Formal-Linear-Like-Factorization 
of a given polynomial. This we can do through the algorithm we describe below. 
It is based on a methodology that is analogous to the usual Euchdean division 
among polynomials and can be implemented on a computer via proper software 
(MATHEMATICA, for instance). 

THE FORMAL ALGORITHM 

Input: A multivariable polynomial p, a set of undetermined parameters W = {Wi^a^^p}, 

taking values in R. 

Initial Conditions: A; = 

Step 1: We set A; = -I- 1. 

Step 2: We find the maximum term of p, maxterm{p) = CkX''^'^ ■ ■ ■ x''n'''. The coefficient 
Cfc, in the first iteration, is a constant number. Then, it depends on the set of 
parameters W. 

Step 3: We form the linear polynomials: ^ = xi, ^2,^ = W2,i,fea;i + X2 L^^k = 

Wn,l,kXl + Wn,2,kX2 + Wn,3,kX3 H h a^n- 

Step 4: We make the subtraction: Rj. = p - CkL{^'^ L^^'j^ ■ ■ ■ 
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Step 5: IF Rj. depends only on the variable xi THEN set R = and go to the output 
ELSE put p = and go to step 1. 

Output: The quantities R, c^, Li^^, jj,^. /J. = I, ■ ■ ■ ,k, i = 1, . . . ,n 

The following theorem proves the finiteness and the efficiency of the algorithm. 

Theorem 2.2 The Formal Algorithm terminates after a finite number of steps. 
If R, c^, Li^i^, ji^^, /J, — 1, . . . ,k, i — 1, . . . ,n are its outputs, then: 

k 

FormalLF\p] = ^ c.L^; • Lj; • Lj; • • • L^- + R 

Proof: Let p be a multivariable polynomial and z = ^x°^x'^ ■ ■ ■ x'^" its maximum 
term. We follow the Formal Algorithm step by step. When A; = 1, step 2 will 
give ci = 7 and ji^i = ai,...,jn,i = an- Taking into consideration the above 
values and the construction of the linear polynomials Li i, L2,i, . . . , L^^i, step 4 will 
produce a new polynomial which will not contain the term z. Obviously, the 
maximum term of Ri, which is called Zi, will be ordered lower than z, Z\ -< z, with 
respect to the order raised earlier. By induction, we get for the maximum terms 
z >~ zi )^ Z2 )^ z^ )^ ■ ■ ■ . This nest and the construction of the order will finally 
eliminate all but Xi-terms. This fact guarantees the termination of the algorithm. 
Substituting now reversely, we have successively p — CiLf 'l ■ • ■ L^^i + Ri, Ri — 

C2L{]i ■ ■ ■ ]J^;^ + i?2, • • • , Rk-i = CklJ^i • ■ ■ + R- Combining these we get: 

p — Yl^^jL=i '^M-^il/T ■ ■ ■ -^^"m + which is the Formal-Linear- Like- Factorization upon 

request. □ 

The following example exhibits the function of the algorithm. 

Example 2.2 We have the previous polynomial p = 5xi — 7xiX2 + 11x1X3. We 
want to find a Formal- Linear- Like- Factorization of p. In order to clarify our ideas 
we shall follow the Formal- Algorithm in detail. First, maxterm{p) = 11x1X3, 
here ci = 11 and ji,i = l,ji,2 = 0,ji,3 — 1. The linear polynomials Lj^i are 
Li,i = xi, L2,i = VF2,i,iXi + X2, 1/3,1 = W^^i^iXi + W3,i,22;2 + X3. Then we get 
Ri =p(a^i,a:2,X3)-llLji-L^i-L^i = 5xi-7xi.X2-llxi(VF3jjXi + M'^3_i,2a;2+a;3) = 
5xi — IIVF3, 1,1X1 + (—7 — 1114^3, i,2)xiX2. The new maxterm is (—7 — llW3,i,2)xiX2, 
with C2 = —7 — llW3,i,2, ii,2 = 1; ^2,2 = 1- Thc polynomials Li2 now become 
Li,2 = xi, L2,2 = W2,i,2a;i +X2 and R2 is, i?2 = -Ri - (-7 - llW^3,i,2)a;i(W2,i,2a:i + 
X2) = 5xi -|- (7^2,1,2 — llW^3,i,i + llW^2,i,2W^3,2,i)2;?. R2 Contains only x\-terms 
and the algorithm terminates. Therefore, the Formal Linear-Like Factorization 
of p is p = 11x1(1^3,1,1X1 + W3^i^2X2 + X3)+ (-7 - llVl/3,i,2)a^i(W^2,i,2a;i + 3^2) + 

+ 5X1 + (7W^2,l,2 - 11H^3,1,1 + llW^2,l,2W^3,2,l)^? 

We can take different expressions of a concrete polynomial p, by giving certain 
values to the parameters Wj,^.,^. Such procedures are called evaluations of the 
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FormalLF[p\. The most rigorous approach is the following: Let W = {lVj,o-,,^} be 
the set of the variables that appear in the Formal - Linear - Like - Factorization 
of a given polynomial p. By arranging the parameters in an increasing order we 
form the vector W = (lVi^,CTfe,vpJfc=i,2,.,n- Let r = (0^)^=1,2,. ,n be a vector of real 
numbers which has the same length as the vector W. We say that the parameters 
W follow the rules r and we write W ^ r if the following substitution is valid: 
^«fe,crfe,¥'fc = dk, k = 1,2, . . . ,n. Let M a set of rules, M = {ri, r2, . . . , r^} then 

A k 

FormalLF[p] | = |J {J^ ^fc^t/T " " " + R with W ^ G M} 

i/=i iJ.=i 

The set of substitutions M, may be finite or infinite. 

Example 2.3 We shall work with the polynomial p that appears in the example 
\2.'J\i . Let M = {(—2,1,-1)}. We determine the following order among the 
parameters (W^3,i,i, W^3,2,i, 1^2,1,2)- Then ^3,1,1 = -2,^3,2,1 = l,W^2,i,2 = -1 and 

FormalLF[p] = llxi(— 2xi + X2 + x^) — 18xi(— Xi + X2) + 5xi + Axf. If 
M = {{ip,(f,9)},ip,9 G R then we set W^^i^i = v^, ^3^2,1 = V^, ^^2,1,2 = d and 

FormalLF[p] |^ = llxi{(pxi + (px2 + X3)+ {—7 — ll(p)xi{9xi + X2)+ 5xi + {79 — 
Hip + 119(f)x1. This is just a re-parameterization of the FormalLF[p]. 

In this section we shall examine the positivity of multivariable polynomials via 
the Formal-Linear-Like-Factorizations developed above. This is a well-known is- 
sue and there are and other computational approaches - recently, for instance, in 
pUj - that are mainly based on algorithms that construct sums of squares. Our 
contribution essentially relies on the existence of parameters that provide us with 
a more flexible tool that is also suitable for solving the positivity problem, among 
others. A given polynomial p G R[a;i, X2, ■ ■ ■ , is called positive, if p > 0, for 
every (xi , X2 , • • • , Xn ) G R". Let S be the set of values of the parameters Wi^a;<p, 
such that all but the square terms of p are eliminated and the coefficients of the 
remaining terms are positive real numbers. Then, if iS 7^ the polynomial is posi- 
tive and Formal LF[p] |^ ia a class of sum of squares of p. Here we need to make 
clear that all the calculations are symbolic and not numerical. This permits us to 
study not only specific polynomials with numeric coefficients, but also classes of 
polynomials with parametric coefficients. In any case, the construction of the set 
S can be carried out via the following algorithm: 

THE POS- ALGORITHM 

Input: A multivariable polynomial p, a set of undetermined parameters W = {Wi^a^ip}, 
taking values in R. 

Initial Conditions: 5 = {}, E = {}. 

Step 1: By means of the FORMAL ALGORITHM we get the quantities R, c^, Lj,^, jj,^, 
^ = I, . . . ,k, i = 1, . . . ,n. We denote the coefficients of the remainder R by c^, 
^ = k + 1, . . . ,k + h, h the number of terms of R. 
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step 2: REPEAT FOR n = 1, . . . ,k + h 

IF some of the exponents jj,^, i = 1, . . . , n are odd numbers THEN 

= 0U {c^} ELSE E = EU {c^} 

NEXT II 

Step 3: Find the values of the parameters Wij^^ that eliminate the " odd " coefficients 
and make the "even" coefficients positive. In other words, we construct the set 

S = {Wij^(f, = Xi,j,(f, G R : = and Cp > with C/j, e O and Cp G E}. 

Output: The set S. 

The proof of the following theorem is straightforward. 

Theorem 2.3 Let p G R[xi, X2, ■ ■ ■ , Xn] be a given multivariable polynomial. If S 
is non-void, then p is positive definite, and FormalLF\p] |^ is a family of " sum 
of squares " expressions of p. 

We would like to make the following remarks in connection with the above algo- 
rithm: 

Remark 2.1 (i) We can modify the algorithm so that the whole procedure is ex- 
ecuted " together " with the Formal Algorithm, as opposed to after it. This will 
allow for a quicker implementation of the method, (ii) Let us suppose that we have 
a Formal - Linear - Like - Factorization of a given polynomial p. The coefficients 
of the first terms with odd exponents contain a small number of parameters ( usu- 
ally one or two). This means that they can be eliminated easily for some particular 
values of the W -parameters. By substituting those values into the other terms we 
decrease the number of parameters, thus simplifying the whole computational pro- 
cedure significantly. Clearly, the method is not computationally complex and can 
be carried out normally. Certain algorithms can be used towards this direction, for 
instance the Stetter algorithn, to mention but a few. (Hi) If the output of the POS- 
AIGORITHM is S = ^, this does not mean that the polynomial p is not positive 
or that another sum of squares does not exist. Nevertheless, our approach can be 
extended by using nonlinear " factors ", for instance: Ln,k — Wn,i,kXi + VKn,2,fea^2+ 

\- Wn,{i,2),kXiX2+ • • • + Wn,{n-i,n-i),kXn-i + x^- this casc wc Can obtain a 

further sums of squares that can successfully deal with cases in which the current 
method fails. This will be undertaken in a future study. 

Example 2.4 Let us consider the polynomial p — x^ — 2xy + 6y^ — Ayz + 3z^. Its 

Formal-Linear-Like-Factorization is: 3(2;-|-a;M^3^i^i-|-yM^3^2,i)^ +(2/+'^W^2,i.2)(~6VF3_2,i — 
4) iz+xW3,i,2+yW3,2,2) +a;(6M/3,2,lW^2,l,2+41^2.1,2-6W3,l.l) iz+xW3,i,s+yWs,2,3)+ 

{y + xW2,l,4? (-3W|2,l+6W3,2,2W^3,2,l+4W^3,2,2 + 6)+ x{y + xW2,l,n) (61^2,1,4^^1,2,1- 
6^3,1,1^^3,2,1+ 6VF3,i,2'W3,2,l+6W2,l,2VF3,2,2VF3,2,l- 12^2,1,4^^3,2,2^^3,2,1 -614^2,1,2^^3,2,3^^3,2,1 
-12^^2,1,4 + 4M/3,1,2+ 4^^2,1,2^^3,2,2 -8^^2,1,4^^3,2,2 -4^^2,1,21^3,2,3 +6H^3,l,lW'3,2,3 - 
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2) + R (We do not include the entire remainder because of its size). In order to 
eliminate the first non-square term to appear, we set ^3^2,1 = ~|- This trans- 
forms the factorization as follows: 3 (— y + z -\- xW^^i^iY —6x1^3^1^1(2; + xW^^i^^ + 
yW^,2,3) +f (2/ + xV^^2,i,4)' -lx{y + xW2,i,5) (14Vr2,i,4-6iy3,i,i-9VF3,i,iiy3,2,3 + 3) 
-\-R'. The values ^3,1,1 = and VF2,i,4 = eliminate the other non-square 

terms and the factorization becomes p = 3(— y + z^ +y(~if + v)'^ +Yi^^- 
Thus, the set S upon request is S = {W3_2,i = — |,VF3,i_i = 0, ^2,1,4 = —jg 
all the other parameters Wi^j^^ are free} . The above " sum of squares " expres- 
sion of p guarantees that p is positive. 

3 The Feedback Asymptotic Stabilization 

We are now in a position to apply the entire concept that was raised in earher sec- 
tions to the problem of feedback asymptotically stabilizing a nonlinear system at a 
given equilibrium point. This is a well-known topic that has been discussed exten- 
sively in the literature [Sj, [Ej. In this paper we adopt a computational approach. 
Specifically, let us have the continuous nonlinear system: x = $(x, u), where 
X = {xi, X2, ■ ■ ■ , XnY is the state vector, u = (mi, W2, • • • 1 the input vector and 
$ = [$i(x, u), . . . , $„(x, u)]"^, where <l>j(x, u), z = 1, . . . , n are multivariable poly- 
nomials of ( ) and {ui, . . . ,Um) without free terms. Obviously (x°, u°) = 
(0, 0) is an equilibrium point. Let us denote by {A, B) the linearization pair of this 
nonlinear system "around" the origin, (i.e. A = ^(0,0), i? = ^(0,0)). As is al- 
ready known, if the linear system {A, B)is asymptotically controllable (which 
we can check by testing the sign of the real parts of the uncontrollable eigenvalues. 
If all of these are negative, then the system is asymptotically controllable), then 
the corresponding nonlinear system is locally asymptotically stable at the origin. 
This means that we can find a matrix F such that the feedback-law constructs a 
system x = $(x, F(x)) which is locally asymptotically stable at the point (0,0). 
The main concern of this paper is to calculate this quantity F computationally, 
but with the following alterations: 

1) F does not need to be linear (a matrix) but may also be nonlinear (a poly- 
nomial function). Actually, the problem under examination is that of finding 
a state feedback of the form u = a(x), with a(x) = [ai(x), . . . , am(x)]-^ and 
aj(x), i = 1, . . . ,m multivariable polynomials too, where the corresponding closed- 
loop system x = $(x, a(x))) has a global asymptotically stable equilibrium at 
(xO,uO) = (0,0). 

2) Some of the uncontrollable eigenvalues of {A, B) can have zero real parts. 

3) The stability is not local but global. 

At this point we introduce the algorithm in order to deal with the feedback asymp- 
totic stabilization problem. 

THE FEEDBACK-GAS ALGORITHM 
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Input: The polynomials $i(x, u), a polynomial Lyapunov function L, the degree of the 
feedback law upon request. 



Step 1: We define the feedback law u = a(x), a(x) = [ai(x), a2(x), . . . , ari(x)], with 

n {n,n) 

n=i (n,i2)=(i,i) 



(n,n,...,n) 

(jl,i2,.--,*n)'^*l'^*2 

(il,i2,.--,«fc) = {l,lv,l) 



and a\^^ ■ ■ N unknown parameters taking values in R. 
Step 2: We define the quantity: 

dL 

V = - — ($i(x,a(x)),$2(x,... ,$„(x,a(x)) 

Step 3: By means of the POS-ALGORITHM we construct the set V, consisting of those 
values of the parameters for which V is positive. 

Output: The set V. 

Theorem 3.1 Let us have the nonlinear continuous system: x = $(x, u). Let V 
he the output of the FEEDBACK-GAS- Algorithm. IfV 7^ then the set of feedback 
laws u = a(x) |^ make the origin globally asymptotically stable. 

Proof: The positive definiteness of the quantity V guarantees that the Lyapunov 
function L decreases along the trajectories of the closed-loop system #(x, a(x)) 
and, therefore, the origin is asymptotically stable. Given that this is the case for ev- 
ery point, the origin is globally asymptotically stable. 



Example 3.1 The angular momentum of a rigid body controlled by two indepen- 
dent torques can be described, following some simplification, through the equa- 

$i(x,u 

with 



tions: x = $(x, u) 



and $1 



$2( 
$3( 



U 



[ui{t),U2{t)] 

03X1X2. The quantities 01,02 



X, uj 
X, u) ^ 

O1X2X3 + Ui, <^2 = a2XiX3 + U2, $ 

and 03 are certain constants and 03 7^ 0, (nl4-h P^9^ 176). It can be proved 
that this system can be globally stabilized about x = and u = 0. A spe- 
cific feedback can be constructed by following certain methods I14]l - Following the 
steps of the Feedback- GAS- Algorithm we select a Lyapunov function of the form 
L = x\ + x\ + x\, and a pair of feedback-laws of the forms: ui = AiXi -\- B1X2 + 
T1X3 + A1X2X3, U2 = A2X1 -\- B2X2 + T2X3 + A2X1X3. Then, we define the quantity: 
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V = -^i{Xi,X2,X3,AiXi + BiX2 + TiX:i+AiX2X3)Xi -^2{Xl, X2, X3, A2X1 + B2X2 + 

r2a:3 + A2a:2X3)x2 —$3(0:1, X2, X3)x3. The Formal Algorithm will give the following 
Formal- Linear- Like- Factorization of V : 

(ai + 02 + as + Al + A2)[VF2,i,i(x3 + XiW^^i^s + X2W3^2,3)xl 

+ {X2 + XiW2,lj){W3^1,l + W2,l,lW3,2,l " 21^2,1,5^^3,2,1 " ^^2,1,1^^3,2,3)2:? + 
{X2 + Xi 14^2,1,5)^(^3,2,12^1 - X2- Xi 1^2,1,1) (2:3 + 2^1^3,1,1 + 2:2^3,2,1)2:1] + 
(12^2,1,2 - ri)(x3 + Xi 1^3,1,4 + 2:2^/3,2,4)2:1 + (^2 + XiVF2,l,8) 
{-A2 -B,+ 2S2V^2,1,6 + r2W^3,l,2 + r2Vl/2,l,2H^3,2,2 " 2r2iy2,l,6W^3,2,2 + TiW3,2,4- 



r2W^2,l,2H/3,2,4)2:i + ix2 + 2;iK/2,l,6) (r2H/3,2,2 " 
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r2(a;2 + 2:1^/2,1,2) (2:3 + 2:1^3,1,2 + 2:2^3,2,2) + R 

where R is the remainder. The values of the parameters that eliminate the non- 
squares terms of the above expressions are M/2,1,1 = 0, VF3,2,i = 0, W3,i,i = 0, 
H/2,1,4 = Ti = 0,T2 = 0, Ai = -A2-ai-a2-a3 and Ai, A2, Si, S2, A2, ar- 

bitrary. For those specific values the FormalLF[V] becomes: ("^2+2-Bi^2+gi 4:AiB2)x^ 

— B2 {^2 + ^^ijf^ ^i) ■ -^'^ order to get positive coefficients we further demand 

B2 < and Ai < ^'^^'^^^^^^^^ ■ These conditions mean that the Lyapunov func- 
tion decreases along the trajectories of the system, which guarantees the asymp- 
totic stability of the origin. The family of feedback laws is given by the rela- 
tions ui = AiXi + B1X2 + AiXi2;2, U2 = A2X1 + i?22;2 + A2X1X3 with B2 < 0, 
Ai < ^^"'"^fjj^^"^^^ ; Ai = A2 — di — a2 — 03 and A2,Bi, A2 arbitrary. We can also 
obtain other classes of feedback, by choosing other expressions for the quantities 
Ui, U2 or other Lyapunov function. 

Example 3.2 Let us suppose that we have the nonlinear system: x = $(x, u) = 



$i(x,u) 
$2(x,u) 



with $1 = Ax + 8y- ^x'^ + 5yui-^UiU2 - Au2 +^xu2, $2 = - ^2: - 

?>y^ — \u\x, where x = [x{t),y{t)Y is the state of the system, consisting of two 
functions and u = [ui{t).,U2{t)\ the input vector. We can easily check that the ori- 
gin is an unstable equilibrium point for the system. Furthermore, the linearization 
of this system is not asymptotically controllable, since the polynomial Xu has a root 
with zero real part (see fT^ for details). Following the steps of the Feedback-GAS 
Algorithm, we choose a pair of feedback laws of the form: ui = Aix + Biy + Tixy, 
U2 = A2X + B2y + T2xy. Then, we define the quantity V = — $i(x, y, Aix + Biy + 
Tixy, A2X + B2y + T2xy)x -$2(2:, y, Aix + Biy + Tixy, A2X + B2y + T2xy)y. The 
Formal Algorithm will give the following Formal- Linear- Like- Factorization ofV: 

FormalLF[V] = -^TlWl^^.x' - hlW2,i,2 (^y + ^^^^2,1,2) x'+ 
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-4r2A2 - 2orir2A2 - 25r2r2 , , , i . 



12^2 ' ' ' ' 2 

+^(-^2 + 32^2^2 + 8oriB2 - leoFi + 8or2Bi)- 

Bl- 576AI - 2880A2ri - 2880^112 + 576012 \ ^ 2 
^ ~ ^72(-S| + 32^252 + 80rifi2 - leOFi + 80r25i) J 

( (-2A2-5ri)\^ ,5, 
+r2U-x^ Ib^ ) B2X^ + 2^B,B2-2B,y 

(-20^2^1 + 40Ai - 3212 - 20^2^1 + 40^2 - 55) V 
^""^ 40B^{B2-2) ; 

/ ^2 \4 

+4(S2-2)(y + a;iy2,i,i3)a; + 3 fx^ + yj +R 

where R is the remainder (not included because of its size). The values of the 
parameters that eliminate the non square terms of the above expression and make 
the coefficients of the remaining terms positive, are V = {Ai = 2, A2 = 

B2 = 2, T2 = 0, with < Bl < I and — q^q^^ ^ — - — - < Ti < 

i6B^ 225Bi+32ViTT-^5B^ B^ y Wjj^o- arbitrary. By taking, for instance, Bi — 
IjFi — J the FormalLF[V] becomes: 

I 20375 4 / 659 V o o (x \4 
Farnu.lLF[V\ |, = — + 2 . + , j + + 3 (- + 

which means that the Lyapunov function decreases along the trajectories of the 
system. This guarantees the asymptotic stability of the origin. The family of the 
feedback laws is given by the relations Ui = x + Biy + Tixy, U2 = -^x + 2y, with 



16Bj -225L'i -:-!2 vTlT ,/r,i,';> -/ji^' ^ F ^ -225i;i -:-!2v'TTTy/o«j^-«« 



u < -Di < 4, ana ^^^r^ < i i < ^^^^ 



It should be noted that the values and | for the parameter Bi are a kind of 
bifurcation point for the feedback expression, above or below which the asymptotic 
stabilization is not achieved. 
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